Bose-Einstein Condensation and Spin Mixtures of Optically Trapped Metastable 
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We report the realization of a BEC of metastable helium-4 atoms ( 4 He*) in an all optical potential. 
Up to 10 s spin polarized 4 He* atoms are condensed in an optical dipole trap formed from a single, 
focused, vertically propagating far off-resonance laser beam. The vertical trap geometry is chosen to 
best match the resolution characteristics of a delay-line anode micro-channel plate detector capable 
of registering single He* atoms. We also confirm the instability of certain spin state combinations 
of 4 He* to two-body inelastic processes, which necessarily affects the scope of future experiments 
using optically trapped spin mixtures. In order to better quantify this constraint, we measure spin 
state resolved two-body inelastic loss rate coefficients in the optical trap. 

PACS numbers: 37.10.Gh, 03.75.Hh, 34.50.-s, 34.50.Fa 



When a helium atom is placed in the long lived 2 3 Si 
metastable state, it behaves in many ways like any num- 
ber of other atomic species used in trapping and cooling 
experiments. In particular, metastable helium (He*) can 
be manipulated by common laser and evaporative cooling 
techniques, the application of the which have culminated 
in the Bose-Einstein condensation of 4 He* [iHj] and the 
realization of a degenerate Fermi gas of 3 He* |5[. Despite 
many similarities however, ultracold metastable helium 
also presents unique rewards and challenges in relation 
to the more traditional ultracold alkali atomic systems 
due to the almost 20 eV of potential energy per atom 
that is stored in the metastable excited state. 

One of the most appealing aspects of He* is the range 
of possibilities it offers for the detection of individual 
atoms, as well as for the real-time monitoring of the 
trapped gas. Through the use of a micro-channel plate 
(MCP) and delay-line anode, it is possible to measure 3- 
dimensional (3D) time of flight (TOF) distributions with 
single atom resolution Q . In addition, the stream of ions 
produced by the collision induced ionization of trapped 
atoms provides a realtime in-situ readout of the density 
of the trapped gas. Using these techniques, we have pre- 
viously studied phenomena such as the onset of a BEC in 
real-time 0] , the Hanbury Brown- Twiss effect for bosonic 
and fermionic atoms 0, and the degenerate 4- wave 
mixing of matter waves resulting from collisions between 
condensates 0, [l(| • The pairing of an MCP and a phos- 
phor screen has also allowed the direct imaging of He* 
atom laser modes (ll- 13|. 



The large energy stored in the metastable state that 
enables the exceptional detection schemes used in these 
experiments comes at a price, however: the loss of atoms 
in the trapped gas throu gh t he enhancement of collisional 
ionization processes 

0,1 Il5j . Indeed, were it not for the 
fact that such losses are strongly suppressed in a polar- 



ized gas EMI, it would not have been possible to reach 
the atomic densities sufficient for condensation. 

To date, experiments involving degenerate He* have 
been largely based upon magnetic trapping, if for no 
other reason than the magnetic trap guarantees the spin 
polarization (for the case of 4 He*) that enables the high 
densities that are necessary to reach degeneracy. In re- 
cent years however, the power and flexibility of all optical 
potentials formed from far-detuned laser beams, together 
with the lifting of the constraint of spin polarization, has 
facilitated access to diverse fields of study including mod- 
els of condensed matter systems and spinor physics [lij . 
The combination of custom engineered trapping poten- 
tials, the spin degree of freedom, and the power of single 
atom detection of He* should open many exciting possi- 
bilities. 

As a first step in this direction, we have constructed 
an optical dipole trap for He*. In this paper, we describe 
the initial implementation of this trap and the subsequent 
production of a BEC. Moreover, we report measurements 
of the dominant ionizing loss rate constants between dif- 
ferent spin states, the values of which will determine the 
feasibility of future studies involving non-polarized gases. 



BEC OF 4 HE* IN AN OPTICAL TRAP 

Up to N = 10 8 4 He* atoms at a temperature of T - 300 
/iK are initially loaded from a magneto-optical trap into 
a Ioffc-Pritchard type magnetic trap in the 2 3 Si, mj = 1 
magnetic substate where they are cooled to T = 150 fjK 
by a retro-reflected ID Doppler beam [l9[ . An RF forced 
evaporation ramp is then applied for ~ 4 s to further cool 
the gas to T = 15 /iK, N ~ 5 • 10 6 , after which time the 
optical trap, formed from a single far detuned 1.5 /im 
wavelength laser beam focussed to a 1 /e 2 gaussian beam 
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waist radius of w — 43 /Ltm, is ramped up to a full power 
of P ~ 1.5 W over approximately 1 s. The magnetic 
trap is then switched off in a time on the order of 100 
/is, and a bias magnetic field B — 4 G is applied in the 
horizontal direction in order to maintain the polarization 
of the atoms. 

At full power, the trap provides U ~ 50 fiK of confine- 
ment in the radial (horizontal) direction, though due to 
gravity, an effective "lip" is formed at the bottom of the 
axial (vertical) potential. The height of this lip and cor- 
respondingly the effective trap depth, U e ff, is strongly 
dependent on the optical trap laser power as well as on 
residual vertical magnetic field gradients. Under typical 
conditions and at full optical power, the trap depth is 
reduced in this direction by a factor of three. To further 
characterize the trapping potential, trap oscillation fre- 
quencies are measured via the resonant loss of the atoms 
resulting from parametric excitations driven by intensity 
modulation of the laser. At full trap depth, the approxi- 
mated harmonic trap oscillation frequencies are o; r /27r = 
2.4 kHz and lo z /2tt = 15 Hz. At this trap depth, we also 
measure a background pressure limited 1 /e atom lifetime 
of t ~ 25 s using a dilute polarized mj = 1 gas, for which 
two and three-body losses are highly suppressed. 

After transfer from the magnetic trap, approximately 
N = 10 6 atoms remain in the optical trap at a temper- 
ature of T = 5 /jK. After 1 s of thermalization, N is re- 
duced by 10% and the temperature reaches a steady state 
value of 3 fiK. Further cooling is achieved by evaporation 
as the trap laser intensity, /, is decreased from maximum, 
I 0l towards a constant non-zero value, If, in an exponen- 
tial fashion over ~ 4 s as I(t) = (I a — I F )e^ t / T " > + I F , 
where the time constant r = 1 s. After evaporation, the 
trapping beam is switched off and the atoms fall onto the 
MCP detector, located 47 cm below the trapping region 
(TOF ~ 300 ms), as illustrated in fig. 0] 

The final degeneracy of the atomic cloud can be con- 
trolled by either varying the offset of the evaporation tra- 
jectory, Ip, or by interrupting the trajectory before the 
final trap depth is reached. Figure [2] shows time of flight 
momentum distributions corresponding to three different 
final trap depths. The distribution of the condensed frac- 
tion exemplifies the striking inversion of aspect ratio of 
the initially cigar shaped condensate that results from 
the anisotropic expansion during the long TOF. 

Due to the vertical orientation of the beam relative to 
gravity, the escaping atoms exit the trap primarily over 
the lip at the bottom of the trap, along the propagation 
direction of the laser, during evaporation. This results in 
an intense beam of optically guided atoms that is regis- 
tered on the MCP directly below the trap. This signal 
of evaporated atom flux vs. time, shown in fig. [2] has 
proven useful in the optimization of atom transfer from 
the magnetic trap and subsequent evaporation, and pro- 
vides a unique realtime method of monitoring the evolu- 
tion of the cooling gas. Moreover, the same mechanism 




FIG. 1: Trap and detector layout. The optical trap is 
formed from a vertically propagating red-detuned beam. The 
anisotropy of the potential results in a vertically elongated 
condensate with typical radial and axial (vertical) Thomas- 
Fermi radii of R r ~ 2.5 fim, and R z ~ 600 /xm, where R\ 
= 2/i3_o / / mijJi (the 3D chemical potential fX3n is defined in 
the text, and trap frequencies u)i corresponding to a BEC are 
given in the following figure caption). After release from the 
trap, during the 300 ms TOF, the condensate expands into a 
horizontally oriented "pancake" distribution since the radial 
size increases by a factor of ~ 3000, while the axial size in- 
creases by only a few percent. After their TOF, individual 
atoms are detected in the horizontal plane with a crossed de- 
lay line anode positioned beneath an MCP. Combined with 
the time of arrival, this allows a 3D reconstruction of the TOF 
atomic distribution. 



has recently been used to guide an atom laser formed 
from a He* BEC [l3|. 



The primary motivation for the vertical orientation of 
this trap, however, is to create a better "mode matching" 
between the time of flight density distribution of the con- 
densate and the resolution capabilities of the MCP de- 
tector, for which the resolution in the vertical direction, 
a z fa 3 nm, is much higher than that of the horizon- 
tal direction, a x , y ~ 300/im. In addition to aligning the 
small dimension of the "pancake" BEC with the high 
resolution direction of the detector (see fig. [1]), this ge- 
ometry is expected to improve the contrast of measure- 
ments of Hanbury Brown-Twiss type correlations, since 
the shortest correlation length, associated with the verti- 
cally elongated potential, is now aligned with the vertical 
axis 

Finally, we note that condensates formed at the low- 
est trap depths, such as those shown in fig. [21 approach 
the crossover to ID that occurs when the 3D chemical 
potential falls below the energy of the first harmonic 
oscillator state of the radial potential [2(J 21 1: — 
\Tiuj{lbNa/a) 2 / 5 < huj r , where Q = (uj^uJz) 1 / 3 is the geo- 
metric mean of the trap frequencies, a — y/h/(mui) is the 
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leads to increased susceptibility to inelastic loss processes 
whose rates are enhanced by Penning and associative ion- 
ization: 





FIG. 2: Evaporation to BEC. On the left is the flux of atoms 
exiting the trap over the vertical potential lip during evapo- 
ration. These atoms remain confined by the trapping beam 
in the horizontal plane, and so are guided onto the MCP. 
The peak near time = s is associated with atoms lost dur- 
ing transfer from the magnetic trap and the subsequent rapid 
thermalization of the remaining atoms. The peak near 4 s is 
the TOF density distribution of the cold atoms released at 
the end of evaporation. In the center are 2D column densities 
obtained by integration along one horizontal axis of the 3D 
TOF density distributions of released clouds of various tem- 
peratures. On the right are the corresponding axial profiles 
made by integration along the remaining horizontal axis. At 
full trap depth, the temperature T — 3 fiK > T c , where the 
critical temperature for condensation T c w 2 fiK, and so the 
momentum distribution is thermal (top). As the tempera- 
ture is lowered below T c , the distribution becomes a bimodal 
superposition of thermal and condensed atoms, marked by 
the high visibility of the pancake shaped condensed fraction 
contrasted with the uniformly distributed thermal component 
(middle). When the trap intensity is lowered until / m 0.5 I , 
the temperature T <C T c , and the thermal component is re- 
moved, leaving a quasi-pure condensate of up to N = 10 5 
atoms with condensed atom fraction N c /N > 0.9 (bottom). 
At this trap depth (U e ff/kB ~ 0.25 /LtK), the trap frequencies 
are LU r /2n — 1.5 kHz and Lj z /2n — 6.5 Hz. 



characteristic oscillator length, a = 7.5 nm is the s-wave 
scattering length [22j, and m is the 4 He atomic mass. 
Condensates with further reduced number N ~ 10 4 re- 
main well within our detection capabilities, and can in 
fact cross this boundary. 



INELASTIC LOSSES IN A THERMAL SPIN 
MIXTURE OF 4 HE* 



As already mentioned, it is primarily the stored in- 
ternal energy of the metastable 4 He* atom that distin- 
guishes it from more commonly trapped alkali atoms. 
This characteristic allows single atom detection, but also 



He* + He* 



He + He+ 
Ret + e~ 



(1) 



In an external field, the 2 3 Si state of 4 He* has three 
magnetic substates mj = ±1 and mj = 0. In a mag- 
netic trap, the required polarization of the trapped gas 
in the m,/ = 1 state suppresses these processes. In an 
optical trap, however, atoms in all spin projections can 
be confined. To better comprehend the feasibility of ex- 
ploiting this added degree of freedom, we now turn to 
state-resolved measurements of inelastic loss from spin 
mixtures of thermal atoms in the optical trap. 

At low temperatures such as those considered in this 
work, where s-wave atomic collisions dominate, inelastic 
collisions between 4 He* atoms can take place through the 
5 E+ potential, for colliding atoms with total spin M ^ 0, 
or the 1 S+ potential for M = 0, where M = Em; is the 
sum of the spin projections of the colliding atoms. Ion- 
ization processes are spin-forbidden for collisions through 
the former, which is the origin of the suppression of in- 
elastic losses in spin polarized samples. Moreover, any 
other combination of input states with M/0 should be 
likewise suppressed. The collisions that remain are those 
between atoms with M = 0, namely collisions involving 
two mj = atoms or one each of mj = 1 and mj = -1 

Emm. 

By creating mixtures of various spin states and mon- 
itoring their populations as a function of time, we can 
test this prediction and also quantitatively measure the 
loss rates for the dominant processes. As previously de- 
scribed, atoms are initially loaded into the optical trap 
and cooled in the 2 3 Si, mj = 1 magnetic substate, the 
same as is trapped in the magnetic trap. After transfer, 
the optical trap depth is lowered such that the gas tem- 
perature approaches, but does not go below, T c . This 
cooling phase is followed by an adiabatic recompression 
to full intensity, where the trap is held for ~ 1 s to ensure 
thermalization. Typical values following this trap trajec- 
tory are, N = 5 x 10 5 , and T = 2.0 fxK, so that T/T c > 
1.25 with peak density n ~ 2.5 x 10 12 cm -3 . 

Following cooling and recompression, spin mixtures are 
created from the polarized mj — 1 gas by a non-adiabatic 
RF sweep across the transition between adjacent levels, 
mj =±1h mj = 0. For the current work, we use 2 ms 
linear sweeps of 3 MHz centered at ~ 11 MHz at three 
separate power settings to prepare the atom spin states in 
different proportions. Figure [3] shows a numeric solution 
to a three level Landau-Zcncr (L-Z) problem as a function 
of sweep strength, which is proportional to the RF power 
divided by the linear sweep rate (25[. As illustrated in fig. 
131 for the case of the weak sweep "a" , the final population 
is approximately evenly split between the mj = +1 and 
mj = states, with a small amount (~ 10%) of mj = 
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FIG. 3: RF sweep schemes. The large panel shows the so- 
lution to the three level L-Z system for input state mj — 1, 
plotted vs. "Sweep Strength" which is proportional to the RF 
power for a sufficiently broad sweep of fixed duration. Panels 
[i-v] show the effects of ramps of various strengths (labeled 
"a", "b", and "c") on different input states. Note that the 
total population is an incoherent mixture of spin states, so 
only single state inputs are considered. Panels [i-iii] corre- 
spond directly to the main panel, and [v] is equivalent to [iii] 
with the mj — ±1 states reversed. The behavior of panel 
[iv] can not be inferred from the main panel, since the input 
state in this case is m.j = 0, though we have confirmed the 
result both experimentally and numerically, as described in 
the text. 
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FIG. 4: Loss processes in imbalanced mixtures. Two ini- 
tial spin state populations prepared using "a" and "b" of fig. 
|3] are plotted vs. time to determine the stable and unstable 
spin combinations. The error bars correspond to the standard 
deviation of the mean derived from repeated measurements. 
The solid curve through the mj — points is a fit to a two- 
body loss process corresponding to a preliminary value /3qq ehm 
~ 7 x I0 -10 cm 3 /s, and is shown for reference, since in this 
case, systematic uncertainties in atom number and tempera- 
ture that result from residual distortions associated with the 
Stern-Gerlach field limit the precision of the measurement. 
The dotted horizontal lines are shown for reference. 



-1. A sweep of intermediate power, "c," places half the 
atoms in the mj — state, with the remaining half split 
evenly between the mj — ± 1 states. As the RF power 
is increased further, the sweep "b" splits the population 
between toj = and toj = -1, with a small amount (~ 
10%) of raj = 1 remaining. 

Before making quantitative measurements of the loss 
rates, we first check the overall relative stability of var- 
ious combinations of spins. After spin preparation, the 
trap is held at full power for a variable time, and is then 
suddenly switched off. A Stern-Gerlach vertical magnetic 
field gradient is then applied in order to impart a spin- 
dependent momentum component, such that after TOF, 
the arrival times of the three spin state populations are 
different and each can be counted separately by the MCP. 

Figure |4^i and |4Jd show measurements of losses in im- 
balanced spin mixtures prepared with sweeps "a" and 
"b" of fig. EJ respectively. In fig. H^, the populations of 
all three spin states initially decay until the to,/ = -1 is 
depleted, at which time the mj = 1 population stabilizes 
while the mj = population continues to decay for the 
entire duration. Figure |4|d shows a similar measurement 
in which the roles of the mj — 1 and mj — -1 have been 
reversed. The observed behavior confirms the instability 
of the (mj,m'j) — (1,-1) and (0,0) mixtures, as well 
as the relative stability of the (m,/,TOj) = (0, 1), (0,-1), 



(1,1), and (-1,-1) mixtures. The extent of the stabiliza- 
tion of the to j = 1 in[4j! and the mj = -1 in [4Jd in the 
presence of the remaining mj = gas places a coarse 
upper limit on the loss rate coefficients /3io,-io,ii,-i-i 
< 10~ 13 cm 3 /s, though two-body loss processes below 
this level are indistinguishable from the background one- 
body losses under the present experimental conditions. 
This finding is consistent with the suppression of ioniz- 
ing collisions already observed for spin polarized mj = 1 
samples, and moreover, is a direct demonstration of the 
extension of this suppression to other atom pairs with 
total spin M ^ 0. 

With the stable and unstable combinations of spin 
states identified, we now focus on a careful and quantita- 
tive treatment of the more dominant loss rates between 
states (mj,m'j) = (0,0) and (1,-1). In the following, 
we model the local atomic density of the thermal gas, 
n = n(r), as 

^ = -rn-(3n 2 -L 3 n\ (2) 

where T = 1/r = 0.04 s^ 1 is the measured one-body loss 
rate due to background gas collisions and off-resonance 
scattering of trap laser photons, j3 is the two-body loss 
rate coefficient, and L3 accounts for three-body contri- 
butions to the loss rate. Though we keep the three-body 
term L3 for the time being, we intentionally use low den- 
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sity gases in order to minimize its contribution and isolate 
two-body effects. 

In order to compare with our data, eqn. [2] is put into 
terms of the total number, N, by spatial integration over 
the extent of the inhomogeneous density distribution of 
the trapped cloud, n(r). 
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^ = -TN - /3N(n) - L 3 N(n 2 ) 
at 



(3) 



where (ra«) = {1/N) J d 3 r[n(r)]« +1 . 

If we consider populations of interacting magnetic 
substates {mj,m'j) = with numbers Ni(j) in the 

i th {j th ) substate, eqn. [JJbecomes 



dN, 



U) 



dt 



-TN 



L 3 N tot {n 2 ot ), (4) 



where j3ij is the partial two-body rate coefficient for loss 
resulting from collisions between states (mj,m' 7 ) = {i,j), 
and (?%•(,)) = (l/ATj^-)) J d 3 r[ni{r)rij (r)] is determined as 
in eqn. [3J by spatial integration of the density distribu- 
tions of both states i, j. Here we have also generalized the 
three-body term as an effective rate constant dependent 
upon the sum of all spin component densities, ntot- 

Figure [5^i shows the loss of raj — atoms from a mix- 
ture prepared with half the atoms in the raj — state 
and the remaining half split evenly between the raj = ± 
1 states as in fig. [3J panel (iii). In order to further sup- 
press the contribution of three-body processes in these 
measurements, we intentionally reduce the atom number 
such that the peak density is reduced to no ~ 10 11 cm" 3 . 
Furthermore, since we have previously shown the stabil- 
ity of the mj = state to loss processes involving the raj 
= ±1 states, we attribute all two-body losses in the raj 
= to interactions among mj — atoms. To improve 
signal to noise on the number and temperature measure- 
ment of the raj = atoms, the uij = ±1 atoms that 
remain after the optical trap hold time are cleared by a 
strong transverse magnetic field gradient that is applied 
shortly after the optical trap is switched off (the raj = 

atoms are not perturbed by this field gradient due to 
their lack of magnetic moment). 

To measure the losses between raj — 1 and raj = - 

1 atoms, we follow a similar procedure with one added 
step. Figure [SJd shows the spin mixture that is initially 
prepared in the same way as in fig. [5^, using sweep "c" of 
fig. [H After the optical trap hold time, however, the RF 
pulse is applied a second time at the same strength be- 
fore the strong clearing magnetic field gradient is applied. 
We find that the second application of the RF sweep at 
power "c" splits the initial raj = population into equal 
parts raj — ± 1, as in fig. [3j panel (iv), which are sub- 
sequently cleared by the field gradient pulse. Likewise, 
the initial raj = ± 1 population is transferred to mj — 
with 50% efficiency (fig. [3] panels (iii), (v)), which pro- 
tects it from the field gradient pulse and other stray field 
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FIG. 5: Loss of mj = and mj = ±1 atoms. The population 
of mj = or the mj = ±1 atoms are separately measured to 
refine measurements of the respective loss rate constants j3oo 
and /3±i. Error bars, corresponding to the standard deviation 
of the mean, are equivalent to or smaller than the data point 
symbols (on the linear scale). The curves show the best fit 
to a two-body decay process for a thermal bose gas in the 
true optical plus gravitational potential, including the mea- 
sured one-body loss rate, as described in the text. The insets 
show the same data and fits plotted as 1/N vs. time. When 
plotted in this way, the decay resulting from purely two-body 
processes appears linear with time, and so the linearity of the 
data exemplifies the limited contribution of one and three- 
body loss processes. 



gradients so that it can fall to the detector unperturbed. 
We have confirmed the observed behavior of the popula- 
tion transfers by further analysis of the three level L-Z 
problem (2(|. Since the populations of raj — 1 and raj 
— -1 are equal, the resulting measured number N±\ — 
\{N 1 +N_ l )=N l =N_ 1 . 

To find the loss rate constants, we numerically solve 
the loss rate eqn. [3J and find the best fit to the data. For 
a thermal gas in a harmonic trap, the atomic density dis- 
tribution is gaussian, and can be analytically integrated 
to give (n) = N/{2y^V) and {a 2 ) = iV 2 /(3v%V 2 ), 
where V — {u> z / u> n)[rauj 2 /2nkBT] 3 / 2 is the volume of 
the trapped thermal gas. 

As expected for our intentionally low peak atomic den- 
sities, no contribution of the three-body term is detected 
in the data. An upper limit L3 < 5 x 10~ 22 cm 6 /s is 
determined by increasing L3 until its effect on the fit to 
the data can not be compensated, to first order, by vary- 
ing the two-body term within its uncertainty. However, 
even before this limit is reached, the quality of the fit 
is significantly degraded. For reference, a value of L3 
«lx I0 -26 cm 6 /s has been measured in a spin polarized 
raj = 1 condensate of 4 He* in a magnetic trap 0, 0, [13] ■ 

By varying the two-body rate constants we find pre- 
liminary values of /3q = 7.6(4) x ICP 10 cm 3 /s and j3 9 ±l = 
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as the num- 

PI- 



8.4(10) x 10~ 10 cm 3 /s, where we define /3q 
ber loss rate coefficient for mj = atoms, (3 g ±l — ^ x _ x 
P 9 _n as the number loss rate coefficient for either mj — 
1 or mj = -1 atoms, and where the superscript g denotes 
the use of the gaussian density approximation. For the 
time being, the uncertainty reflects only the statistical 
and fitting uncertainty. 

We refine the fit values of f3 by integrating eqn. [3] over 
the atomic density distribution of a bose gas at the mea- 
sured temperature in the true optical plus gravitational 
potential. We find that the inclusion of the slight trun- 
cation of the density distribution by the finite trap depth 
[iH and additional terms of the thermal bosonic density 
distribution result in a slight reduction of the final ex- 
tracted values of /3 o = 6.6(4) x 1CT 10 cm 3 /s and /3±i = 
7.4(10) x 10~ 10 cm 3 /s. 

In addition to the statistical errors given above, several 
sources of systematic error must also be included. The 
dominant contributions come from the determination of 
the absolute number of trapped atoms and the character- 
ization of the trapping potential. Together with residual 
uncertainties in temperature [33j these combine to give 
an overall estimate of jjSp sys rj 25%. 

For a comparison of relative rates, the systematic er- 
rors are largely reduced. We find that within the remain- 
ing uncertainty that /3oo = P±i, meaning that all three 
spin components deplete with the same rate constant, 
which in turn confirms that the collision "event rate" for 
the (m'j,mj) — (0,0) collisions is half that of the (1,-1), 
since each (0,0) collision results in the loss of two mj — 



atoms 

Previous measurements of the inelastic loss in unpolar- 
ized He* gases have been performed by measuring ioniza- 
tion rates of samples prepared in magneto-optical traps 
171 l29l432j , with the results of, and discrepancies among, 



these various measurements summarized in refs. |24l32|. 
These experiments differ from that of the present work 
in that the losses resulting from interactions between in- 
dividual spin states are not resolved, and instead an ion- 
ization rate, Ki" npol \ is measured for the global decay of 
the entire unpolarized gas. To compare with these mea- 
surements, our partial two-body loss coefficients can be 
re-scaled to give the global loss rate using 

Pmeas = lSK^^ ~ (2P-1P+1 + pi) (5) 

where (3 mea s is the measured number loss rate, and pi 



is the fractional population of the mj = i substate 17 1 



The numerical factors are defined such that the ioniza- 



tion rate constant K^ npo ^ — ^/3 meas (i.e. 1 ion detected 
per 2 lost atoms, see eq. [1} for the case of an unpolar- 
ized gas where p_i — p+i — Po = ■§• Our measured /3oo 
with p Q = 1 (p ±1 = 0) gives K{T pol) = y 00 = 1-10(7) 



x lO" 10 cm 3 /s. 

3 2 



Similarly, if p-i = p+\ = 



K 



(unpol) 



one half in (3±i reflects the fact that we measure only 
N\ or N-i to determine /3±i, whereas when both states 
are counted, the rate constant is reduced by half. As 
before, a systematic uncertainty of 25% applies to these 
numbers. These values of K^ np are in good agreement 



with those measured in refs. 13 1[ l32l a nd are consistent 
with theoretical predictions 15, 



In conclusion, we have achieved Bose-Einstein conden- 
sation of 4 He* in an optical dipole trap. We observe ef- 
ficient evaporative cooling and a long trap lifetime, both 
of which point favorably towards the possibility of fu- 
ture all-optical production of these condensates as well 
as the implementation of more complex trapping geome- 
tries. Using this trap, we have also made the first spin 
state-resolved measurements of two-body loss in a ther- 
mal 4 He* gas, and find that it is dominated by collisions 
between (mj,m'j) — (0,0) and (1,-1) atoms. We find loss 
rate coefficients for these individual processes that are 
in line with theory and compatible with previous mea- 
surements of the global loss rate from magneto-optical 
traps. Finally, we conjecture that although these rela- 
tively high loss rates may prevent the cooling of some 
non-polarized mixtures to degeneracy, spin mixtures can 
be created starting from a cold polarized gas whose life- 
time will not preclude all experiments. This is especially 
true for the case of highly polarized spin mixtures con- 
sisting of a small number of mj = atoms immersed in 
a gas of mj — 1 or -1 atoms. 
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of the successful operation of an optical dipole trap for 
He* at the VU University in Amsterdam 34 1 



1.23(17) x 10~ 10 cm 3 /s, where the factor of 



[1] A. Robert, O. Sirjean, A. Browaeys, J. Poupard, 
S. Nowak, D. Boiron, C. I. Westbrook, and A. Aspect, 
Science 292, 461 (2001). 

[2] F. Pereira Dos Santos, J. Leonard, J. Wang, C. J. Bar- 
relet, F. Perales, E. Rasel, C. S. Unnikrishnan, M. Leduc, 
and C. Cohen-Tannoudji, Phys. Rev. Lett. 86, 3459 
(2001). 

[3] A. S. Tychkov, T. Jeltes, J. M. McNamara, P. J. J. Tol, 
N. Herschbach, W. Hogervorst, and W. Vassen, Phys. 
Rev. A 73, 031603(R) (2006). 

[4] R. Dall and A. Truscott, Optics Communications 270, 
255 (2007). 

[5] J. M. McNamara, T. Jeltes, A. S. Tychkov, W. Hoger- 
vorst, and W. Vassen, Phys. Rev. Lett. 97, 080404 
(2006). 

[6] M. Schellekens, R. Hoppeler, A. Perrin, J. V. Gomes, 
D. Boiron, A. Aspect, and C. I. Westbrook, Science 310, 
648 (2005). 

[7] S. Seidelin, O. Sirjean, J. V. Gomes, D. Boiron, C. I. 



7 



Westbrook, and A. Aspect, Journal of Optics B: Quan- 
tum and Semiclassical Optics 5, S112 (2003). 
[8] T. Jeltes, J. M. McNamara, W. Hogervorst, W. Vassen, 
V. Krachmalnicoff, M. Schellekens, A. Perrin, H. Chang, 
D. Boiron, A. Aspect, et al., Nature 445, 402 (2007). 
[9] A. Perrin, H. Chang, V. Krachmalnicoff, M. Schellekens, 
D. Boiron, A. Aspect, and C. I. Westbrook, Phys. Rev. 
Lett. 99, 150405 (2007). 

[10] V. Krachmalnicoff, J.-C. Jaskula, M. Bonneau, G. B. 
Partridge, D. Boiron, C. I. Westbrook, P. Deuar, P. Zin, 
M. Trippenbach, and K. Kheruntsyan, (to be published) 
arXiv.org:0911.4564 (2009). 

[11] R. G. Dall, L. J. Byron, A. G. Truscott, G. R. Dennis, 
M. T. Johnsson, M. Jeppesen, and J. J. Hope, Opt. Ex- 
press 15, 17673 (2007). 

[12] R. G. Dall, L. J. Byron, A. G. Truscott, G. R. Den- 
nis, M. T. Johnsson, and J. J. Hope, Phys. Rev. A 79, 
011601(R) (2009). 

[13] R. G. Dall, S. S. Hodgman, M. T. Johnsson, K. G. H. 
Baldwin, and A. G. Truscott, Phys. Rev. A 81, 
011602(R) (2010). 

[14] M. W. Miiller, W. Bussert, M. W. Ruf, H. Hotop, and 
W. Meyer, Phys. Rev. Lett. 59, 2279 (1987). 

[15] P. S. Julienne and F. H. Mies, J. Opt. Soc. Am. B 6, 
2257 (1989). 

[16] G. V. Shlyapnikov, J. T. M. Walraven, U. M. Rahmanov, 
and M. W. Reynolds, Phys. Rev. Lett. 73, 3247 (1994). 

[17] N. Herschbach, P. J. J. Tol, W. Hogervorst, and 
W. Vassen, Phys. Rev. A 61, 050702(R) (2000). 

[18] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008). 

[19] P. O. Schmidt, S. Hensler, J. Werner, T. Binhammer, 
A. Gorlitz, and T. Pfau, J. Opt. Soc. Am. B 20, 960 
(2003). 

[20] F. Schreck, L. Khaykovich, K. L. Corwin, G. Ferrari, 
T. Bourdel, J. Cubizolles, and C. Salomon, Phys. Rev. 



Lett. 87, 080403 (2001). 
[21] A. Gorlitz, J. M. Vogels, A. E. Leanhardt, C. Raman, 

T. L. Gustavson, J. R. Abo-Shaeer, A. P. Chikkatur, 

S. Gupta, S. Inouye, T. Rosenband, et al., Phys. Rev. 

Lett. 87, 130402 (2001). 
[22] S. Moal, M. Portier, J. Kim, J. Dugue, U. D. Rapol, 

M. Leduc, and C. Cohen- Tannoudji, Phys. Rev. Lett. 

96, 023203 (2006). 
[23] V. Vcnturi and I. B. Whittingham, Phys. Rev. A 61, 

060703(R) (2000). 
[24] P. J. Leo, V. Venturi, I. B. Whittingham, and J. F. Babb, 

Phys. Rev. A 64, 042710 (2001). 
[25] M.-O. Mewes, M. R. Andrews, D. M. Kurn, D. S. Durfee, 

C. G. Townsend, and W. Ketterle, Phys. Rev. Lett. 78, 

582 (1997). 

[26] G. E. Carroll and F. T. Hioe, J. Opt. Soc. Am. B 2, 1355 
(1985). 

[27] O. Sirjean, S. Seidelin, J. V. Gomes, D. Boiron, C. I. 

Westbrook, A. Aspect, and G. V. Shlyapnikov, Phys. 

Rev. Lett. 89, 220406 (2002). 
[28] O. J. Luiten, M. W. Reynolds, and J. T. M. Walraven, 

Phys. Rev. A 53, 381 (1996). 
[29] H. C. Mastwijk, J. W. Thomsen, P. van der Straten, and 

A. Niehaus, Phys. Rev. Lett. 80, 5516 (1998). 
[30] M. Kumakura and N. Morita, Phys. Rev. Lett. 82, 2848 

(1999). 

[31] P. J. J. Tol, N. Herschbach, E. A. Hessels, W. Hogervorst, 
and W. Vassen, Phys. Rev. A 60, R761 (1999). 

[32] R. J. W. Stas, J. M. McNamara, W. Hogervorst, and 
W. Vassen, Phys. Rev. A 73, 032713 (2006). 

[33] We compensate, to first order, for a small (~ 10%) sys- 
tematic increase in temperature that occurs as the gas 
decays by rescaling the measured number to account for 
the resultant slight decrease in peak atomic density. 

[34] W. Vassen, private communication. 



